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A uniform dilute Bose gas of known density has a universal behavior as the atomic scattering 
length tends to infinity at unitarity while most of its properties are determined by a universal 
parameter ^ relating the energies of the noninteracting and unitary gases. The usual mean-field 
equation is not valid in this limit and beyond mean-field corrections become important. We use 
a dynamical model including such corrections to investigate a trapped disk-shaped dipolar Bose- 
Einstein condensate (BEC) and a dipolar BEC vortex for large scattering length. We study the 
sensitivity of our results on the parameter ^ and discuss the possibility of extracting the value of 
this parameter from experimental observables. 
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I. INTRODUCTION 

The properties of a uniform dilute interacting — Bose 
or Fermi — atomic quantum gas interacting by an S'-wave 
contact inteaction at zero temperature is determined by 
two scales — the atomic scattering length a and den- 
sity n. As a — > oo at unitarity, the first scale is not of 
concern and the observables of the gas are solely deter- 
mined by density and the gas exhibits a universal be- 
havior. Bulk chemical potential p, of the unitary gas 
is proportional to the only available energy scale — the 
Fermi energy (or the chemical potential of the noninter- 
acting gas) Ep — h'^{6'K^n)^^^ /2m ~ so that p. — ^Ep, 
where ^ is a universal parameter and m is the mass of 
an atom [l|-[l6j. Similarly, the energy per particle E of 
the unitary gas is proportional to the energy per particle 
E{= 3Ef/5) of the noninteracting gas: E = ^E 0, g. 
The Fermi energy is a physically meaningful quantity for 
the Fermi gas, but the same can also be used as an en- 
ergy scale for the Bose gas 0, [l^, [II]. The Bose and 
Fermi gases behave similarly at unitarity, because the 
Bose gas exhibits fermionization. If this fermionization 
of the Bose gas is absolute, then ^ should be the same 
for the Bose and Fermi gases. 

The parameter ^ relating the energy of the noninter- 
acting and unitary gases, has been "measured" experi- 
mentally from a study at or near unitarity of the density 
1 , ' O'" of ground-state energy, or of sound velocity 
8| of a trapped Fermi gas. But a similar experiment 
is more difficult for a Bose-Einstein condensate (BEC) 
due to a large probability of three-body loss by molecule 
formation at or near unitarity, which is the threshold 
for molecule formation |12| . In the weak-coupling limit 



(ani/3 < 1) the mean-field Gross-Pitaevskii (GP) equa- 
tion gives a good description of a trapped BEC, where n 
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is the density. In the strong-coupling regime {an^^^ > 1) 
the GP equation highly overestimates the atomic contact 
interaction and leads to unphysical results. Experimen- 
tal activities to access the strong-coupling regime of a 
BEC, to test the beyond mean-field corrections pTi. [Tsj. 
and to extract the parameter ^ from these studies have 
just began [l2[. Although unitarity is also the thresh- 
old for molecule formation in a two-component (spin up 
and down) Fermi gas, the probability of formation of di- 
atomic molecules is highly suppressed in this case due 
to Pauli repulsion among spin-parallel fermions in the 
three- fermion system and hence is not of concern l[ • 

Lately, BECs of ^^Cr [H and ^^-^Dy [lU, [13 with 
large dipolar interaction have been observed and studied. 
The inter-atomic interaction now has two components: 
an S'-wave contact interaction and an anisotropic long- 
range dipolar interaction. This allows to sttidy the dipo- 
lar BEC with a variable contact interaction [l^ [22| using 
a Feshbach resonance [2^. The intrinsically anisotropic 

tolar BEC [H has many distinct features [H [13, 
. The stability of a dipolar BEC depends not only 
on the scattering length, but also on the trap geome- 
try [l^, [2^ [23]. A disk-shaped trap leads to a repul- 
sive dipolar interaction and the dipolar BEC is more 
stable, whereas a cigar-shaped trap yields an attractive 
dipolar interaction and hence favors a collapse instability 

We study the static and dynamic properties of a disk- 
shaped dipolar BEC and dipolar BEC vortex, with the 
dipole moments aligned perpendicular to the plane of the 
disk, for lar ge s cattering length using a beyond-mean- 
field model |16l . [30j for the BEC-unitarity crossover. 
In this paper we consider the strong-coupling limit of 
the contact interaction only and not the same limit of 
dipolar interaction. In the weak-coupling limit, this 
crossover model reduces to the GP equation and the 
Lee-Huang- Yang (LHY) correction [sij, whereas in the 
strong-coupling regime of large scattering length it re- 
duces to the universal result at unitarity. We find that 
the radial densities are sensitive to the parameter ^ in 
the strong-coupling regime and hence a study of density 
is expected to yield information about this parameter. 
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However, the frequency of oscillation of the dipolar BEC 
is found to be insensitive to this parameter. From a study 
of vortices in a disk-shaped dipolar BEC we find that 
both density and radius of vortex core are sensitive to 
this parameter in the strong-coupling regime. 

In the disk configuration, the dipolar interaction is 
highly repulsive, as parallel dipoles arranged in a plane 
with the dipole moment perpendicular to the plane repeal 
each other [l^, [2^ . The strongly repulsive dipolar inter- 
action should reduce three-body loss by molecule forma- 
tion in the strong-coupling regime, as the rate of the re- 
action 3A^A2+A should be suppressed in this setting 
with A representing a dipolar atom and A2 a molecule. 
Also, as both the contact and long-range dipolar inter- 
actions contribute to molecule formation, the threshold 
for molecule formation will be displaced from unitarity, 
specially in strongly dipolar BECs, thus creating a new 
scenario of experiment with a dipolar BEC in the strong- 
coupling regime to determine the parameter ^. 

In Sec. II we present the mean-field and beyond- 
mean-field models to study a dipolar BEC in the weak- 
and strong-coupling regimes as well as along the BEC- 
unitarity crossover as the scattering length a is increased. 
We also present a Gaussian variational formulation for its 
solution at unitarity. In Sec. Ill we present the results of 
numerical and variational studies of density, root-mean 
square (rms) sizes, chemical potential, and frequencies of 
radial and axial oscillations of a disk-shaped dipolar BEC 
and BEC vortex. Finally in Sec. IV we present a brief 
summary and conclusion. 



II. ANALYTICAL CONSIDERATION 

A. Dipolar Gross-Pitaevskii Equation 

We consider a disk-shaped dipolar BEC of N atoms, 
each of mass m, using the CP equation [l^ 



— + V + fl{a,N)+F{add,N) 



with the bulk chemical potential 

/2(a, N) = Airan, 



(1) 
(2) 



and density n = A^|(?!)p. Here the dipolar nonlinearity 
F{add,N)^Nj C/rfd(r-r')|0(r',t)|2dr', 



V 



(3) 



is the harmonic trap, r — {x,y,z} = {p,z}, 
UddCR-) = 3add{l - 3cos^e)/R^, R = r - r', normal- 
ization J 0(r)^c?r ~ 1, 9 the angle between R and the 
polarization direction z. A 3> 1 the trap anisotropy, 
add = fJ.ofj,^m/{12'!Th^) the strength of dipolar interaction, 
fiQ the permeability of free space, and p. the (magnetic) 



dipole moment. In Eq. ([T|) the r and t dependence of p, 
and F arc not explicitly shown and length is measured in 
units of Iq = y^h/mu!, where cu is the harmonic trap fre- 
quency in a; or y directions, time t in units of to = uj~^. 
At unitarity, the bulk chemical potential of Eq. ([T]) is 
independent of a and is 



(4) 



To obtain a quantized vortex of unit angular momen- 
tum h] around z axis, we introduce a phase (equal to the 
azimuthal angle) in the wave function [33 |. This proce- 
dure introduces a centrifugal term l/[2(a-^ + y^)] in the 
potential of the CP equation so that 



V=lix'+y' + X'z')+ / 



(5) 



We adopt this procedure to study an axially-symmetric 
vortex in a disk-shaped dipolar BEC. 



B. BEC-Unitarity Crossover 

Lee, Huang, and Yang (LHY) [3l| obtained the leading 
terms of the beyond-mean-field expression for energy of 
a uniform Bose gas from which the following ex pre ssion 
for the bulk chemical potential can be obtained |33j : 



/i(a, N) = Anan 1 + aa^''^y/n 



(6) 



The lowest order term in this expansion is the CP result 
^ first derived by Lenz [slj. However, expression 
although gives the leading correction for larger an^^^, 
diverges in the strong-coupling regime, and hence has 
only limited validity along the BEC-unitarity crossover. 

In addition to studying the system in the weak- 
coupling limit ^ and unitarity dH), we also consider 
the system along the full BEC-unitarity crossover from 
weak to strong coupling, as the parameter an^^^ is in- 
creased. For this purpose we consider the following mini- 
mal crossover model for the bulk chemical potential con 
sistent with weak and strong couplings 



;ai pote 



^-i(a, TV) = 47rn2/3/(x), x = an 

fix)- ^ + (1 + ^)"^^^^ 



1/3 



(7) 
(8) 



where a = 32/(30?), a/7 = ^(67r2)2/3/87r, and v is the 
only free parameter in this expression. The parameters a 
and 7 are determined by the constraints that expression 
([7]) be consistent with the LHY correction ^ as well 
as the unitarity limit both independent of the pa- 
rameter ^. Expression ([7|) is weakly sensitive to v and 
a smooth interpolation between the weak and strong- 
coupling regimes is obtained for any small v. In this 
study we use v ~ 1. This value of v was used Q suc- 
cessfully in a study of ^Li2 BEC in the BEC-unitarity 
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crossover. Expression ([7]) can also reproduce fairly well 
[30| the energies of diffusion Monte Carlo (DMC) calcu- 
lation [sH of a trapped bosonic system of small number 
of atoms. It was found that the energies obtained from 
the crossover model (O for that bosonic system are in 
better agreement with the DMC calculation than those 
obtained from the GP equation. 

The crossover model ([7]) is Galilei invariant and yields 
the hydrodynamic equations of the dipolar BEC at zero 
temperature, and enables one to study collective dynam- 
ical properties of the system in the full crossover from 
weak-coupling to unitarity [l^ [s^]- Equations ([T]) and 
(O should be considered as a generalization of the GP 
equation with beyond mean-field corrections to properly 
include the effect of interaction for large positive scatter- 
ing length. The saturation of the interaction at unitarity 
is properly taken care of in the crossover model ([7]). As 
an application we shall study here the properties of a 
disk-shaped dipolar BEC and dipolar BEC vortex in the 
strong-coupling regime to show the sensitivity of the re- 
sult to the universal parameter ^. 



effective Lagrangian L = J Cdr becomes 



L = - [2wla 



wl6 



III 



27r ifpWz 



9C 



2/3 4/3 ' 
Wz Wp 



(11) 



where C = ^/i^{Q■n'^Nf/'^/(b{)■n^/b), k = Wp/w^, and 



1 + 2^2 SK^tanh" Vl - 



1 



1 



(12) 



The corresponding Euler-Lagrange equations govern- 
ing the evolution of the widths Wp and Wz yield 



add Ng{K) 
%/27r wlwz 



UC 



2/3 7/3' 
Wz Wp' 



... ,2 1 "dd 2iV 

Wz + \ Wz = ^ 7= o y h{K) 



12C 



2tt w^w^ 



5/3 4/3- 
Wz Wp 



where 



(13a) 
(13b) 



C. Variational Approximation at Unitarity 

At unitarity, the dipolar mean- field equations ([l}, ([3|), 
and (U]) can be conveniently solved by a time-dependent 
Lagrangian variational approach. This can be used to 
study the size and frequencies of oscillation of the dipo- 
lar BEC at unitarity. This is done by reducing Eq. ([1]) 
to a system of second order nonlinear ordinary differen- 
tial equations involving the variational parameters. The 
Lagrangian density of Eq. ([T]) is given by [28| 



0>t) + ^|V</.p + ^(p^ + A2z^)|0|2 



Udd{r-v'Mr')fdr'. 

(9) 



5(«) = 



/i(k) 



2 - 7k2 - gK^tanh" Vl - 



1-H10k2-2k"^ QK^tanh^Vl - 



(1 - K 



2~i2 



(1-^2). 



(14a) 
(14b) 



Equations (|13ap and (jl3bp provide the dynamics of the 
evolution of radial and axial widths, respectively. One 
can obtain the expression for the frequencies and lowest- 
lying modes from these equations [38|- The widths for a 
stationary state can be obtained by setting Wp — and 
iuz = in Eqs. (|13a|) and (jl3bp . The chemical potential 
fi for the stationary state is given by 

2\Wp 2w^ J i ^ ^ 
^Nadd /(At) , 1 



27r wf,Wz 



5 9C 

q 2/3 4/3 ■ 
-J Wz Wp 



(15) 



Recalling that n = N\4>\'^, it can be straightforwardly 
verified that Eqs. ([Ij and (|4|) are the Euler-Lagrange 
equations for the Lagrangian density ^ [3p, which 
should be used in the variational formulation [37|. To 
develop the variational approximation, we consider the 
following Gaussian ansatz for the wave function j28j 



(r,t) = 



W^Wz 



exp 



2wl 



z 

2^ 



tap 



(10) 



where the time-dependent variational parameters Wp and 
Wz are the radial and axial widths and a and S are the 
chirps. The Lagrangian density can be calculated by sub- 
stituting the wave function ([TUl) in Eq. Then the 



III. NUMERICAL CALCULATION 

We perform numerical simulation of the 3D GP equa- 
tion ([1]) using the split-step Crank- Nicolson method [33] ■ 
The evaluation of the dipolar integral term in this equa- 
tion in coordinate space is not straightforward due to the 
divergence at short distances. However, this has been 
tackled by evaluating the dipolar term in the momentum 
(k) space. The integral can be simplified in Fourier space 
by means of convolution as [25| 



dr'Uddir - r')|0(r')|' = J-H^[M(k)^[|0n(k)}(r), 

(16) 

where ] and ^{} are the Fourier transform (FT) and 
inverse FT, respectively. The FT of the dipole potential 
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is known analytically [25|. The FT of density is eval- 
uated numerically by means of a standard fast FT (FFT) 
algorithm. The dipolar integral in Eq. ([T]) involving the 
FT of density multiplied by FT of dipolar interaction is 
evaluated by the convolution theorem (|16p . The inverse 
FT is taken by means of a standard FFT algorithm. The 
FFT algorithm is carried out in Cartesian coordinates 
and hence the GP equation is solved in three dimensions 
irrespective of the symmetry of the trapping potential. 
In the Crank-Nicolson algorithm we used space step 0.1, 
time step 0.002 and employed upto 512 space discretiza- 
tion points in each Cartesian direction. We made an 
error analysis of the results for chemical potential and 
rms sizes and found that the maximum numerical error 
of the results reported here is less than 0.5 %. 



A. Experimental Considerations 

Of the experimental dipolar BECs — ^^Cr and ^^*Dy — 
realized so far, the magnetic moment of ^^Cr is /i = 6/iB 
p^ . where ub is the Bohr magneton, and that of ^^^Dy is 
/i = lO^B [24I. Consequently, add = /xo/2^m/(127rfi.^) = 
15ao for ^^Cr and add = 130ao for ^^^Dy, with oq the 
Bohr radius. Hence the dipolar interaction in ^^'*Dy is 
about 9 times stronger than in ^^Cr and we employ a 
is^Dy BEC in this study. For ^^"^Xy^tm estimate for 
the scattering length is a ~ IOOoq [24I. In the actual 
experiment on ^^''^Dy a dipolar BEC of 15000 atoms in 
a fully anisotropic trap with frequencies {fx,fy,fz} = 
{380,500,1570} Hz was obtained [H]. In this study, 
to simulate this experiment [2^ . we use the frequencies 
{fxjyjz} = {436,436,1570} Hz, so that A = 3.601, 
where we take a geometrical mean of the frequencies in x 
and y directions to generate an axially-symmetric BEC. 
The le ngth sca le employed here, for ui = 2tt x 436 Hz, is 
lo = ^/h/Jmuj) = 0.376 yum. 

Next we summarize the different theoretical and ex- 
perimental estimates of ^ obtained so far for bosons and 
fermions. Often the parameter ^ is written as ^ = {1 + (3) 
and different estimates of /3 is given in Table U where the 
variational calculations of Refs. [l^ [l5j for bosons are 
upper bounds and the experimental result of Ref. [l^ 
for ^Li is a lower bound only. Yet another estimate of 
^ can be obtained from a consideration of Fermi super- 
fluid in the BEC side of the Bardeen-Cooper-Schrieffer- 
BEC (BCS-BEC) crossover. Here we reconsider an anal- 
ysis [§] of the experiment on ^Li in the BEC side 
of BEC-unitarity crossover. The molecular BEC of ^Li2 
was then studied using Eqs. ((1]), ([7]), and (|H]) but with 
a/7 = Cmoi(67r^)^/V(27r), where ^moi is the universal pa- 
rameter of Ref. Hi, in place of a/7 = 5(67r^)^/^/(87r) 
considered here. This implies that for a comparison of 
the two studies we should take ^ = 4^nioi- The analysis 
of Ref. [1| yielded ^moi ~ 0.4, so that ^ « 1.6 correspond- 
ing to /3 ss 0.6. In that analysis it was assumed that 
the bosonic molecular unitarity of ^Li2 was achieved for 
the same strength of atomic interaction as the fermionic 



TABLE I: Theoretical result and experimental evaluation of 
the parameter /3 = — 1) for the Bose and Fermi gases. 
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unitarity of ^Li. Actually, the bosonic molecular unitar- 
ity should be achieved at a different value of interaction 
and into the BEC side of the BCS-BEC crossover, where 
the system is less repulsive. This would lead to a smaller 
value of the parameter ^ < 1.6) and f3. Hence the 
analysis of Ref. Q gives an upper bound. From the re- 
sults reported in Table |ll the most accurate theoretical 
3 and experimental [1, [l^ estimates for a Fermi gas 
converge to a value of ^ very close to 0.4. 



B. Disk-shaped dipolar Bose-Einstein condensate 

We study a disk-shaped dipolar BEC of 15000 ^^'^Dy 
atoms with add — 130ao as in the experiment of Lu et 
al. [l^. The parameter ^ can be extracted from the 
observables of the dipolar BEC in the strong-coupling 
regime where the observables would be sensitive to this 
parameter. For this purpose, in this paper, in addition to 
the numerical study at unitarity, we also present a com- 
plete numerical study of the dipolar BEC in the strong- 
coupling regime for scattering length a > lOOao using the 
crossover model ^ . 

In Fig. [T] (a), we plot the radial density of the 
BEC along x axis \^{x,y = 0)p obtained by inte- 
grating out the z dependence of density: \^{x,y)\'^ = 
J dz\(j>{x,y, z)\'^ . In this figure we show the result for 
a = lOOflQ using the GP model ^ and at unitarity 
(|4]) in addition to the results for a ~ 100ao,300ao, and 
500ao using the BEC-unitarity crossover model ^ with 
^ = 0.4 and ly = 1. For small a, the density from the 
crossover model ([7]) is in agreement with the GP model 
([2]) and hence practically independent of the parameter 
^, whereas for large a it approximates the unitarity limit 
(III) with the increase of a. In Fig. [1] (b) we plot the radial 
density for a = 500ao for different ^ using the crossover 
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FIG. 1: (Color online) (a) Radial density along x axis 
\^(x,y = 0)1^ = Jdz\(l){x,y = 0, z)f of a "''Dy EEC with 
N = 15000, A = 3.6, add = 130ao,iy = 1 and ^ = 0.4 for 
a — lOOao, 200ao, lOOOao and at unitarity using the GP limit 
(l2|), crossover model and at unitarity ((4]). (b) The same 
for A = 10, a = 500ao and ^ = 0.2, 0.4, 0.6 using the crossover 
model (III). The density at unitarity Q for ^ — 0.4 is also 
shown, (c) The same for A = 10 and different ^ at unitarity 
0. (d) The same for ^ = 0.4, A = 3.6, 10 at unitarity Q for 
add = 130ao (DBEC), and add = (BEG). 



model ([7]). The result at unitarity (jj)) for ^ = 0.4 is also 
shown. The density is sensitive to the parameter ^ for 
a = 500ao as can be seen from Fig. [1] (b) comparing 
the results for ^ = 0.2, 0.4 and 0.6. The sensitivity of the 
density on ^ at unitarity is illustrated in Fig. [T](c), where 
we show the radial density for ^ = 0.3,0.4,0.5 and 1.6. 
Finally, in Fig. [l](d) we show the density at unitarity for 
nondipolar and dipolar BECs for two values of the trap 



FIG. 2: (Color online) (a) The numerical rms sizes (x), (z), 
and chemical potential fi of the dipolar BEG with = 
15000, Odd = 130ao,5 — 0.4,;/ = 1, A = 3.6 versus akp using 
the crossover model (O, as well as at unitarity Q (arrow), 
(b) The numerical Q and variational (v) results for the same 
at unitarity versus add- (c) The numerical Q and variational 
(v) results for the same at unitarity versus (d) The numer- 
ical Q and variational (v) results for the same at unitarity 
versus A. 



asymmetry A = 3.6 and 10. The difference between the 
two densities is more pronounced for A = 10, where the 
dipolar repulsion is stronger. 

In Fig. [2] (a) we plot the rms sizes of the BEC (x), (z), 
and the chemical potential fi for ^ = 0.4, and add = 130ao 
as calculated using crossover model (O, versus the di- 
mensionless parameter kpa, where kp is the Fermi wave 
vector in a harmonic trap defined by kp = (48A^)^/^/r 
[H, where r= y^h/mu!,ui = 27r(/^/y/^)i/3, {f^^jyjz} = 
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{436,436, 1570} Hz. One can find from Fig. [2] (a) liow 
these quantities — (x), (z), ^ — approacfi tlicir values 
at unitarity as a increases. With the increase of a these 
quantities saturate rapidly to their respective values at 
unitarity. In Fig. [2] (b), we plot the numerical (n) and 
variational (v) results for (x), (z), and /i versus add at 
unitarity, which shows that the results are weakly sensi- 
tive to a variation of add- In Fig- ID (c) we plot (x), (z), 
and yLt versus ^ at unitarity, which shows that the results 
are sensitive to a variation of ^. Finally, in Fig. [5] (d) 
we plot (x) , (z) , and fi versus A at unitarity showing the 
strong sensitivity of the results to a variation of A. From 
Figs. [2] (b), (c), (d) we see that the variational results 
are in good agreement with the numerical ones. 

In a recent experiment, Navon et al. (l^ were able to 
make measurements for densities of a very dilute BEG 
of ^Li for a = 2150ao and extract the parameter £ from 
a theoretical analysis using the LHY correction [3l|, . 
The very dilute BEG prepared in a weak trap, allowed to 
make an experiment for large a = 2150ao. But because of 
the low density, the BEG remained away from the strong- 
coupling regime even for a = 2150ao and was studied by 
the LHY correction, rather than a full crossover model 
as in the present study. In this regime the densities are 
weakly sensitive to the parameter ^ and only an upper 
limit ^ < 0.6 could be obtained from that study p^ . 

Apart from density profile and rms sizes of the dipo- 
lar BEG, other observables which can be studied are the 
frequencies of radial and axial oscillations Q,p and O^, 
respectively, of the fundamental modes. We calculated 
these frequencies by numerically solving the variational 
equations (|13ap and (jl3bp in different cases. The ini- 
tial widths Wp and Wz were taken as their equilibrium 
static values and their time evolution is obtained. The 
frequencies of radial and axial oscillations were extracted 
from the time evolution of the respective widths. In 
Fig. [3] (a) we show these frequencies for ^ = 0.4 ver- 
sus add and in Fig. [3] (b) we plot these frequencies for 
add = 130ao versus ^ for A = 3.6, 10. We also calculated 
the axial frequency fi^ from the small oscillation of the 
rms axial size of the BEG upon real time evolution of 
the mean- field equations ([T]) and Because of a mix- 
ture of frequencies of higher modes, it was not possible 
to obtain precisely the frequency fip from a solution of 
the mean-field equations. The mean-field and the vari- 
ational results for Q.^ are in good agreement with each 
other. These frequencies are practically insensitive to a 
variation of add as well as of ^. Hence it may not be very 
fruitful to study these frequencies in the strong-coupling 
regime in order to extract the parameter ^, specially for 
a moderate density as in this study. 
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FIG. 3: (Color online) (a) Variational frequency (v) of radial 
and axial oscillation Sip and fiz of the dipolar BEC with TV = 
15000, ^ = 0.4, V = 1,\ = 3.6, 10 versus add at unitarity from 
Eqs. (|13ap and H13b|) . The numerical frequencies (n) of axial 
oscillation obtained from a solution of Eqs. ([1]) and ((7| are 
also shown, (b) The same for add ~ 130ao versus ^. 



the parameter ^. In this case the radius of the vortex core 
is an observable directly related to the healing length [l[ 
of the BEG and will also be considered. The radius of 
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C. Dipolar BEC Vortex 

Next we study the density of a disk-shaped dipo- 
lar BEG vortex of unit angular momentum for strong- 
coupling and demonstrate the sensitivity of the result on 



FIG. 4: (Color online) (a) Radial density along x axis 
\^ix,y = 0)f of a ^^*Dy BEC vortex of iV = 1000, A = 3.6, 
add = 130ao,i' = 1 and ^ = 0.4 for a = 100ao,200ao, lOOOao 
and at unitarity using the crossover model and at unitar- 
ity 0. (b) The same for ^ = 0.4, A = 3.6, 10 at unitarity 
for add = 130ao (DBEC), and add = (BEC). 
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rg versus afej? for A'' = 1000, Odd = 130ao, A = 3.6,^ = 
0.4, = 1 obtained using the crossover model ([7]). The 
resuh at unitarity Q is also shown. The relative radius 
of vortex core reduces with the increase of the scatter- 
ing length. Similar reduction of the radius of vortex core 
was predicted for a nondipolar BEC before [l^ |4l| . In 
In Fig. E] (b) we plot (z), rp versus ^ at unitarity 

dH) for = 1000, A = 3.6, Odd = 130ao,A = 3.6, = 1. 
As scattering length a increases in Fig. [5] (a) or the pa- 
rameter ^ increases in Fig. [S] (b), the system becomes 
more repulsive leading to a smaller healing length. Con- 
sequently, the relative radius of vortex core ro, which is 
closely related to the healing length, decreases [l|. The 
relative radius of vortex core shows much sensitivity to 
the scattering length a and ^. 



IV. CONCLUSION 



FIG. 5: (Color online) (a) The rms sizes (x), (z), chemical 
potential /i, and relative radius of vortex core ro of the dipolar 
BEC vortex with N = 1000, aad = ISOao,^ = 0.4, u ^ 1, \ ^ 
3.6 versus akp using the crossover model (O, as well as at 
unitarity Q (arrow), (b) The same at unitarity versus ^. 

the vortex core pc is defined as the radial distance from 
the center of the vortex to a point where the density in- 
creases to the maximum value. It is more appropriate 
to consider the relative radius of vortex core defined by 
To = pc/ {x) , which gives the vortex core radius in relation 
to the radial size of the condensate. It is demonstrated 
that the relative vortex core radius rg could be sensitive 
to ^ in the strong-coupling regime and could be useful in 
deciding the value of ^. 

In Fig. m (a) , we plot the radial density of the BEC 
vortex along the x axis |$(a;,y = 0)p for A^ = 1000, 
A = 3.6, Odd = 130ao,^ = 0.4, i/ — 1 for different values 
of scattering length a. In this figure we show the result 
for a ~ lOOflQ and at unitarity (U) in addition to the 
results for a = IOOoq, 200ao, and lOOOao using the BEC- 
unitarity crossover model ([7]). For small a, the density 
obtained using the crossover model ([7]) is in agreement 
with the GP equation ^ and hence independent of the 
parameter ^, whereas for large a it approximates the uni- 
tarity limit In Fig. |l](b), we show the radial density 
at unitarity for nondipolar and dipolar BECs for two val- 
ues of the trap asymmetry A = 3.6 and 10. The difference 
between the two densities is more pronounced for A = 10, 
where the dipolar repulsion is stronger. 

In Fig. [S](a) we plot chemical potential p, and rms sizes 
(x), (z) together with the relative radius of vortex core 



The properties of a BEC at unitarity is controlled by 
a universal parameter ^ relating the energies of nonin- 
teracting and unitary uniform gases. Using the BEC- 
unitarity crossover model ([7|) we studied the properties 
of a disk-shaped dipolar BEC and dipolar BEC vortex 
in the strong-coupling regime. We find that the density 
profiles are sensitive to the parameter ^ in this regime 
and a study of density should yield an information about 
this parameter. We also studied the frequencies of the 
fundamental modes of radial and axial oscillation of this 
BEC and find that they are not much sensitive to ^. For 
a dipolar BEC vortex, in addition to density, the rela- 
tive radius of vortex core is also found to be sensitive to 
^ in the strong-coupling regime, so that a study of this 
radius may reveal information about ^. Also to extract 
the parameter ^ it is not necessary to study the system 
at unitarity. The density profile of the BEC is sensi- 
tive to the parameter ^ for the contact interaction lying 
between the weak-coupling GP and strong-coupling uni- 
tarity limits, so that a study in this regime should reveal 
information about this parameter. In this study we used 
a dipolar BEC of 15000 ^^^Dy atoms in a disk-shaped 
trap of anisotropy A = 3.6, as in the experiment of Ref. 
j22| . and also A = 10. For an experimental study the 
anisotropy of A = 10, or larger, and a BEC with strong 
dipolc interaction is to be preferred. 
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